IDENTITIES OF EULERAIN AND ORDERED STIRLING 
NUMBERS OVER A MULTISET 



JOON YOP LEE 

Abstract. By considering Eulerian numbers and ordered Stirling num- 
bers of the second and third kinds over a multiset, we generalize identi- 
ties of Eulerian numbers and Stirling numbers of the second and third 
kinds and provide g-analogs of these generalizations. Using these gen- 
(jj 1 eralizations, we also compute Eulerian numbers and ordered Stirling 



numbers of the second kind over a multiset. 



1. Introduction 



o 
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Let /^y be an Eulerian number, which is the number of permutations of 
{1, 2,...,d} with i descents. A well-known identity of Eulerian numbers is 
Worpitzky's identity [1] 
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and the following Carlitz's identity [21 [3] is a (/-analog of this identity 
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where [Z]q = IL>=i(l ~ q x ~ m+i )/(l — (f) and A^i(q) is a polynomial of q. 
^ ■ Eulerian numbers are defined over an ordinary set, thus one way to generalize 

identities ([T]) and (J2J) is to consider Eulerian numbers over a multiset. By 
considering Eulerian numbers over a multiset, we may be able to generalize 
identities ©and ©. 

We can apply such an idea to identities of Stirling numbers of the second 
and third kinds. Let {^} be a Stirling number of the second kind, which is 
the number of partitions of a d element set into k nonempty sets. String 
numbers of the second kind satisfy the following polynomial identity 

(3) * d = X>! 
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A Stirling number of the third kind [1] (or a Lah number), denoted by 
, , is the number of ways to partition a d element set into k nonempty 
linearly ordered subsets For example, the partitions of {1,2,3} into 

2 nonempty linearly ordered sets are {(1),(2,3)}, {(1),(3,2)}, {(2), (1,3)}, 
{(2), (3, 1)}, {(3), (1, 2)}, and {(3), (2, 1)}. Lah [7] defined Stirling numbers 
of the third kind by the following polynomial identity 



<4) 4 x -\ +i )^ 

v J k=l 



By considering Stirling numbers of the second and third kinds over a mul- 
tiset, we may also be able to generalize identities ([3]) and and provide 
q- analogs of these generalizations. 

Motivated by these ideas, we will generalize identities (pQ), ©, and dU 
and provide (/-analogs of these generalizations. In particular, we generalize 
identities ([3]) and © and obtain g-analogs of these generalizations by con- 
sidering ordered Stirling numbers of the second kind and third kind 



over a multiset, rather than ordinary Stirling numbers of the second 
and third kinds. Using generalizations of identities ([1]) and ([3]), we also 
compute Eulerian numbers and ordered Stirling numbers of the second kind 
over a multiset. 

The gist of our idea is the following. For a sequence of finite sets of lattice 
points So, S\, 5*2, ... , we compute the numbers of elements in S n by two 
different ways and obtain a polynomial identity. To obtain a q-analog of this 
identity, we compute the following generating function 



I +■ * •■_■-! — I-.'',; 



q 

{x 1 ,x 2 ,...,x d )es 71 



by two different ways. 



2. A TRI ANGULATION OF THE PRODUCT OF SIMPLEXES 

In this section, we introduce a triangulation of the product of simplexes. 
This triangulation will be the main tool to generalize identities (JTJ) , ([3]) , and 
([!]) and provide g-analogs of these generalizations. 
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For a positive integer n, let [n] = {1, 2, . . . , n} and [n]o = {0, 1, . . . , n}. 
We denote a point in M. d by x = (x\, X2, • • • , Xd)- We denote the zero vector 
of M. d by e^o and for each i G [d] we denote the ith unit vector of M, d 
by &d,i = (0, . . . , 0, 1, 0, . . . , 0) (with 1 in the ith coordinate). We define 
ckq to be a d-dimensional simplex whose vertexes are e^o, e^o + e^i, . . . , 
ed,o + e<j,i + ■ ■ ■ + ed,d- Note that «q is the set of points x such that 1 > 

X\ > X2 > ■ ■ ■ > Xd > 0. 

Let d = d\ + d,2 H Vdi be a sum of nonnegative integers. For each j G [I] 

writing a point of M. d i by Xj = (xj : i,Xj : 2, ■ ■ ■ , Xj,dj), we also denote a point of 
R d = n 5=i K^' by x = nj=i xj. Letting {(^,4) | ft G [d]} = {J l j=1 {j} x 
for a point x in R d we always assume that 

1- — x j2,«2 — ' ' " — x jd,idi 

2 - if = x jh+i,i h +i then eitner < 3h+i or = (jh+i,h+i - !)■ 
Let d = (di , c?2 , . . . , ) . We denote the product of simplexes oft , oft , . . . , 
aft by a d = ]J l j=1 oft and the vertexes of a d by e^,...,,, = n^i^o + 

H h e^.,^) where (ii,i 2 , • • • ranges over the set 11^=1 Kj]o- 

Let 5(d) = S(di,d 2 , be a multiset such that for each j G [I] the 

number of j in 5(d) is dj. We define 6(d) to be the permutation set of 5(d) 
and denote a permutation of 5(d) by a = [a±, a 2 , ■ ■ ■ , For a permutation 
a of 5(d) we define o, d {ij ) to be a ci-simplex with the vertexes &i h 1: i h 2 ,--.,ih i 
for ft G [d] such that (i ,i, i ,2, • • • , *<y) = e i>0 and (i fti i, i h>2 , . . . , ih,l) = 
Yla=i e i,<?a f° r h £ [^]- By definition, a d ((j) is the set of points x such that 

L X <Jl,i\ — f (72,42 — ' ' ' — X &d^d1 

2 - if = X ^h+i,ih+i tnen eitner °h < °h+i or {a h ,i h ) = (a h+1 ,i h+ i - 1). 

Let x be a point of a d . Letting Xj 0: i = 1 and Xj d+1 ,i d+l = 0, we can 
represent x by the following sum 

d 

X = ~^2( X jh,ih ~ X jh+l,ih+l) e ih,l,ih,2,---,ih,l 
k=0 

where (io,i,*o,2, • • • ,*o,«) = e i.o and (i M , i h<2 , . . . , i h j) = Yla=i e i,j a for h G 
[d] . This implies that x is a convex sum of the vertexes of a d ( [ j\ , j 2 , . . . , jd] ) • 
Since this is true for every point of a d , it follows that a d is the union of 
the ci-simplexes a d (a) for a G 6(d). Moreover, the set composed of the 



4 JOON YOP LEE 

d-simplexes a d (o~) for a € ©(d) and their faces is a triangulation of a d by 
definition. We denote this triangulation by T a a. 



3. EULERIAN NUMBERS OVER A MULTISET 

Let R be a subset of M. d . For a nonnegative integer n we denote nR = 
{nr | r € R} and we define Z(J2) to be the set of lattice points in R. Recall 
that a lattice point is a point with nonnegative integers entries. For a point 
x let g x = q Xl+X2 ~\ <~x d _ To generalize identity ([I]), we will compute the 
number of elements in Z(na d ) by two different ways and to obtain a q- 
analog of this generalization, which is a generalization of identity ([2|), we 
will compute fi(q) = !C x eZ(na d ) ^ by ^ wo different ways. 

For a permutation cr of S(d) let -D(cr) be the descent set of a, that is, the 
set of indexes h such that > <Jh+l- We define A{a) to be a subset of a d (a) 
composed of points x such that x Uh ^ h > x ah+1 ,i h+1 for all h where h ranges 
over -D(cj). Then for each point x in a d there is a unique permutation a such 
that A(a) contains x. Therefore a d is the disjoint union of the sets A(o~) for 
a € ©(d). We call this disjoint union of a d the first decomposition of a d . 
Stanley [8] obtained such a decomposition when d\ = d% = ■ ■ ■ = di = 1, 
that is, a d is an /-dimensional hypercube. 

The set Z(na d ) is the product of Z(na dl ), Z(na d2 ), . . . , Z(na dl ) and for 
each j G [/] the number of elements in the set Z(na^) is ( n ~^ d] ), thus 

|Z ( na d )| = n|Z(na^)l = n( n ^'). 

A point x is an element of 7 J (nA(a)) if and only x is a lattice point such 
that 

2. 35<r h) i h > x cr h+1 ,4 +1 + 1 whenever h is a descent of <r. 

Thus the number of elements in r L{nA{a)) is ( n_ l- D ^ cr )l+ d ) . Therefore if we 
denote by / d \ the number of permutations of S(d) with z descents, then 
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by the first decomposition of a d we obtain 

\Z(na d )\ = |Z( l+J nA(a))\ = J] \L{nA(a))\ 
o-ee(d) o-ee(d) 

d-l 



E 

CT ee(d) 



n- |D(a)| +d 



E 



d\ fn — i + d 



As a result, we obtain the following identity 



d\ / n — i + d 
i/\ d 



(5) iir d -i=E 

i=i V «W i=0 

Note that Kim and Lee [Polytope numbers and their properties, arXiv:1206.0511 j 
derived identity ([5]) by using the concept of polytope numbers. Since identity 
([5]) is true for every nonnegative integer n, it follows that 

d-i 



nrr j 



E 



d\ ( x — i + d 
d 



j=l y "■> ' i=0 

Theorem 3.1. Eulerain numbers over a multiset S(d) satisfies 



JJ ( x + d 3 



3=1 



d-l 

E 

8=0 



x — i + d 
d 



When d\ = d% = ■ ■ ■ = di = 1, identity ([6]) becomes identity (pQ). Thus 
identity ([6]) is a generalization of Worpitzky's identity pQ . 

Using identity ([6]), we can obtain Eulerian numbers over S(d) in the 
following way. Let Mi be a d x d matrix defined by 



Mi 







G) 

ft 1 ) ffl 








/2<f-l\ /2d-2\ /d 
V d J \ d ) " ' \a 

Then from identity ([5]) we can obtain the following matrix identity 







" ffi ••• 
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Since the inverse of Mi is 



Mf 1 



o 

(-i) x (T) (-i)Vo 1 ) 



by the Gaussian elimination, it follows that 



d 



d + 1\ J-r (h + d. 



E(-i)"(t + ^)n v d , 

/l=0 V 7 7 = 1 V J 



Theorem 3.2. Eulerian numbers over a multiset S(d) are given by 



B-i)"(?:;)np;*) 



Now we consider fi(q). By definition, x,- is in Z(nag j ) if and only if Xj is a 
lattice point such that n > Xj\ > Xj^ > • • • > > 0, thus the generating 
function V d . o x j is a q-binomial coefficient ""t^ [91 [TO] . Since 

Xj eS(na J ) L a i J q 

Z(na d ) = nLi Z(naf ) and g x = J]Li 9 Xj , it follows that 



(7) 



n + d j 



d i J, 



For a permutation cr of 5(d) we define the major index of a to be 
maj(a) = Ylj^Dta) 3- A lattice point x is in Z(nA(a)) if and only if 



x a h ,i h > X o- h+1 ,i h+1 



for h £ [d + 1] \ L>(ct) 



a^.v, > + 1 for h € 



Therefore 



E 



E 



xEZ(nA(cr)) xeZ(nA(cr)) k=l 



n- \D(a)\ +d 
d 



IDENTITIES OF EULERAIN AND ORDERED STIRLING NUMBERS OVER A MULTISET 

Since Z(na d ) is the disjoint union of the sets Z(nA(a)) for a € 6(d), we 
obtain 



Mi)= E E 

o-e6(d) xeZ(nA(cr)) 

~n - \D(a)\ +d 



(8) 



. ^ qmaj(a) 

aee(d) 
d 

E^ife) 



d 



i=l 



n + 1 
d 



where A^ ti (q) = ^ <res(d) q ma 3\ a ) . Therefore combining identities ([7]) and 

D(cr)|=d-i 

dHL we obtain 



(9) 



n 



n + d j 



n + i 
d 



3 iq i= i 

Identity @ is true for every nonnegative integer n, thus by the substitution 
n — > x we obtain the the following (/-analog of identity (0) 
i 



(10) 



n 



x + d j 
d; 



J -l« i=l 

Theorem 3.3. ^4 q-analog of identity is 



x + i 
d 



J q 



n 

3=1 



x + 
d.j 



E a <m(<?) 



d 



J 9 



J J 9 i=l 

Note that identity (|1U|) is a generalization of Carlitz's identity ([2]). 



4. Ordered Stirling numbers of the second kind over a 

multiset 

Similar to the case of Eulearin numbers over a multiset, to generalize 
identity ([3]), we will compute the number of elements in Z(na d ) by two dif- 
ferent ways and to obtain a g-analog of this generalization, we will compute 
Exez(na d ) <? x by two different ways. 

For two vectors v = (yi, V2, ■ ■ ■ , uj) and v' = (v[,v' 2 , ■ ■ ■ we define 
v < v' if Vi < v[ for all % £ [/]. Let 2T(a d ) be the set of simplexes in T a a that 
contains eo^,...^ and edi,d2,—,dr ^ or a ^-dimensional simplex a fe (k € [d]) in 
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I(a d ) let {ei h 1> i h 2 ,...,i ft ; | h G [k]o} be the vertex set of a k where 

(io,i,k,2, ■ ■ -,io,i) = 
( n ) { (ih,i,h,2,- ■ -,ih,i) < («/.+l,i> i W,2r' -,ih+i,l) for h G [k - 1] Q • 
k (h,i,ik,2, • • • ; ijfc,j) = (rfii ^2) • • • , dj)' 

We define I(a fe ) to be the set of convex sums x = Ylh=o c h e i h i,ih 2,-4 i °f 
the vertexes of a fe such that > for h G [& — 1]. This is implies that if 
x is a point of I(a fc ), then there are k — 1 numbers /ii, /12, • • • , in [d] 

such that 



(12) 



X jh 1 ,ih 1 > 



.1 



^ihfc-i+i^hfe-i+i 



u 3dM 



Let x be a point of a d . Suppose that x satisfies condition (|12p . Then a 
simplex a k in X(a d ) with the vertex set {ej a ^ 2 ,...,i a ; I a £ Wo} satisfying 

(io,i,io,2) • • • ,*o,i) = 

{ia,l,ia,2, ■ ■ ■ ,ia,l) = J2b=l e l,jb for a £ [k - 1] 

is the unique simplex such that I(a k ) contains x. Since I(a k ) is a subset of 
a d , it follows that a d is the disjoint union of the sets I(a k ) for a k G Z(a d ). 
We call this disjoint union of a d the second decomposition of a d 

Let a k be a /c-dimensional simplex in X(a d ) with the vertex set {e^ lj j h 2 ,...,i h , | f< 
[k]o} that satisfies condition (jlip . We define the major index of a fc to be 
maj(a k ) = Ylh=i Y!j=i For each /i G [k] if we denote 5 ft = S(i hj i,i hj2 , h,i 
S(ih-i,i, ih-1,2, ■ ■ ■ , ih—i,l)> then (Si, S2, ■ ■ ■ , Sk) is an ordered partition of 
5(d) into k nonempty multisets. Therefore the number of /c-dimensional 
simplexes in X(a d ) is the number of ordered partitions of S(d) into k 
nonempty multisets. We call this number an ordered Stirling number of 



the second kind and denote it by {fc} - Note that if d\ = c?2 
then { d k } = kl{ l k }. 



di = 1, 
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By definition, the number of elements in the set 1i(nl(a k )) is ( n t 1 ), thus 
by the second decomposition of a d it follows that 



"Q ( n + d J 



3=1 



Y[ \Z(na d i)\ = \Z(na d )\ 



3=1 



jZ(nl+J |+J I(a k ))\ 
k=l a k € z{a d i'-> d l) 

£ £ mm(a k ))\ 

k=l a k £l(a d ) 
d 



n + 1 
k 



As a result, by the substitution n — > x we obtain 



(13) 



n 

3=1 



x + d,- 



,1 I ' 5-/ 1 /• , 

a 3 J k=1 l fc J O 



x + 1 
k 



Theorem 4.1. Ordered Stirling numbers of the second kind over S(d) sat- 
isfies 



nc;*)-£{:} cr 



Similar to the case of Eulerian numbers over a multiset, we can obtain 
ordered Stirling numbers of the second kind over a multiset as follows. Let 
Mi be a d x d matrix defined by 



Mo 



(!) o 
(?) ® 








© ® ••• ffl 



Then we can rewrite identity ()13p in the following form 



l /dj 



n;=i 



l fl+dj 



n; =1 



nj. =1 1 i 



(!) o 

(?) a 



{^} 
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Since by the Gaussian elimination the inverse of M2 is 

(-l)°(o) ••• 



M^ 1 



it follows that 



(-i)°ffl 

(-i^-H/x) (-i) d - 2 (/ 2 ) ••• 



{2} -d-d 



/ =1 \ J 



Theorem 4.2. Stirling numbers of the second kind over S(d) are given by 



o 



h=0 



n t 

7=1 v ^ 



For each a k € ^(a d ) the sum of g x for x G Z(nl(a k )) is q ma ^ ak ) 



n+l 
k 



Thus it follows that 



= E E E 

fe=l a fc eX(o d ) xeZ(n;(a fc )) 

= qinaj(a k ) 
k=l a k el(a d ) 



n+l 
A; 



k=l 



n+l 

k 



where B d>k (q) = Ea*ez( a d) <7 m ° i(a&) - In addition, / a (g) = nj=i 
therefore we obtain 

1 



n+dj 



n 

3=1 



n + 



n+l 

k 



1 k=i 

As a result, by the substitution n — > x we get the following <j- analog of 
identity (JT3]) 



(14) 



n 

3=1 



x + dj 
dj 



1 k=i 

Theorem 4.3. A q-analog of identity H3\) is 



I 



n 

3=1 



x + d j 
dj 



E S d,fc(9) 



9 fe=l 



X + 1 



X + 1 

fc 
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Note that identity (|14p is a generalization of the following identity [IT] 



X>^-l) ad)fc(g) 



J 1 k=0 

where a^kio) is a polynomial of q. 



5. Ordered Stirling numbers of the third kind over a multiset 

To derive a multiset version of identity we will compute So-s6(d) \%(na a 
by two different ways and to obtain a g-analog of this identity, we will com- 
pute f 2 (q) = Eaee(d) £xez(na>)) 9 X b y two different ways. 

Let x be a point in Z(na d ). Since a d is the disjoint union of the sets I(a k ) 
for a k € Z(a d ), there is a unique simplex a k in I(a d ) such that Z,(nl(a k )) 
contains x. If a d (o~) is a (i-dimensional simplex in 2T(a d ) that contains a fc , 
then we can represent (a d (o~),a k ) by a /c-tuple of the form 



(15) 



(((7i, . . . ,0-jJ, (cr il+ i, . . .,(7j 2 ), . . . , (cr ifc _ 1 



+1; 



By a simple computation, the number of such fc-tuples is ( di d2 dXfe-l)' 
Note that the vertex set of a fc is {ej h u i h 2j ... t i h l \ h G [k]o} where 

(k,i,k,2, ■ ■ -,io,i) = 

(ih,i,ih,2, • • • , ih,l) = Si=i e;,j for /i G [fc - 1] 
(ik,l,h,2, ■ ■ ■ ,ik,i) = (di,d 2 ,...,di) 

Therefore if we denote by [^Jo the number of /c-tuples of form (|15p . which 
is called an ordered Stirling number of the third kind (or an ordered Lah 
number) over a multiset, then we obtain 

d 

£ |Z(na d (a))| = E E E 



a-ee(d) 



o-e6(d) fe=l o'Ca^ir) 
d 

E 



k=l 



d 




A 





'n + 1 
, k 



In addition, E CT e6(d) l z ( na< V))l = U,^,...,^) )> thus it; follows that 

d 



d \(n + d 
d l ,d 2 , ■ ■ ■ ,dj \ d 



E 

fc=i 



o 



n + 1 

k 
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As a result, we get a multiset version of identity dU) 

d 



(16) 



d \ fx + d 
di,d 2 ,...,dij\ d 



E 

k=l 



O 



x + l 
k 



Theorem 5.1. Ordered Stirling numbers of the third kind over S(d) satisfy 

d 



d \ fx + d 
di,d 2 , . . . ,dj V d 



E 

k=l 



O 



X + l 

k 



Now we consider fzfa). The number of permutations of 5(d) is [ di d ^ d ) 

, thus it follows that 



n+d 
d 



g 




( d ) 


n + d 


\d 1 ,d2,--- ,dij 


d 



and Y;x&(na d (a)) <7 X 

(17) h(q) 
In addition, by the second decomposition of a d we can obtain 

ui)= E E f 

o-ee(d) x£Z(na d (a)) 

= E E E E <? 



crG6(d) fc=l a fe Gl(« d ) x€Z(nI(a fc )) 



E E E 

fc=l<T6S(d) o'6i(o d ) 



maj (ce k ) 



n + 1 

k 



(18) 



E C ^(9) 



fc=i 



n + 1 

k 



where C A>h (q) = Y,ae6(d) Y.^^i^A) q ma ^ a \ Therefore by combining iden- 
tities (fTT|) and (fTHj) we obtain 



d 

d\ , <fo, ■ ■ ■ , di 



n + d 
d 



n + 1 

k 



i k=l 

As a result, we get the following q- analog of identity (|16[) 

d 



d 



d\,d2, ■■■ ,di 



x + d 
d 



i k=i 

Theorem 5.2. A q-analog of identity H16\) is 



x + l 

k 



d 

d\ , d 2 , . . . , d\ 



x + d 
d 



= E C <^) 



x + l 
k 
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